Implementation of a Toffoli Gate with Superconducting Circuits 
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The Toffoli gate is a three-qubit operation that 
inverts the state of a target qubit conditioned on 
the state of two control qubits. It makes universal 
reversible classical computation [lj possible and, 
together with a Hadamard gate |2j, forms a uni- 
versal set of gates in quantum computation. It 
is also a key element in quantum error correc- 
tion schemes [3H3- The Toffoli gate has been 
implemented in nuclear magnetic resonance |3j, 
linear optics |8j and ion trap systems [9|. Ex- 
periments with superconducting qubits have also 
shown significant progress recently: two-qubit al- 
gorithms [10J and two-qubit process tomogra- 
phy have been implemented three-qubit en- 
tangled states have been prepared [12, 13j, first 
steps towards quantum teleportation have been 
taken [14] and work on quantum computing ar- 
chitecture has been done [16j. Implementation of 
the Toffoli gate with only single- and two-qubit 
gates requires six controlled-NOT gates and ten 
single-qubit operations [15j, and has not been re- 
alized in any system owing to current limits on 
coherence. Here we implement a Toffoli gate 
with three superconducting transmon qubits cou- 
pled to a microwave resonator. By exploiting 
the third energy level of the transmon qubits, we 
have significantly reduced the number of elemen- 
tary gates needed for the implementation of the 
Toffoli gate, relative to that required in theoreti- 
cal proposals using only two-level systems. Using 
full process tomography and Monte Carlo pro- 
cess certification, we completely characterized the 
Toffoli gate acting on three independent qubits, 
measuring a fidelity of 68.5 ±0.5 per cent. A simi- 
lar approach [16j realizing characteristic features 
of a Toffoli-class gate has been demonstrated 
with two qubits and a resonator and achieved 
a limited characterization considering only the 
phase fidelity. Our results reinforce the poten- 
tial of macroscopic superconducting qubits for the 
implementation of complex quantum operations 
with the possibility of quantum error correction 
schemes [17J. 

We have implemented a Toffoli gate with three trans- 
mon qubits (A,B and C) dispersively coupled to a mi- 
crowave transmission-line resonator, in a sample which 
is identical to the one used in ref. [T3]. The resonator is 
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FIG. 1. Circuit diagram of the Toffoli gate, a, A not- 
operation (©) is applied to qubit C if the control qubits (A 
and B) are in the ground (o) and excited state (•) respectively. 

b, The Toffoli gate can be decomposed into a CCPHASE gate 
sandwiched between Hadamard gates (H) applied to qubit C. 

c, The CCPHASE gate is implemented as a sequence of a 
qubit-qutrit gate, a two-qubit gate and a second qubit-qutrit 
gate. Each of these gates is realized by tuning the |11) state 
into resonance with 1 20) for a {n, 2n, 3n} coherent rotation 
respectively. For the Toffoli gate, the Hadamard gates are re- 
placed with ±7r/2 rotations about the y axis (represented by 
R^y). d, Pulse sequence used for the implementation of the 
Toffoli gate. During the preparation (I), resonant microwave 
pulses are applied to the qubits on the corresponding gate 
lines. The Toffoli gate (II) is implemented with three flux 
pulses and resonant microwave pulses (colour coded as in c). 
The measurement (III) consists of microwave pulses that turn 
the qubit states to the desired measurement axis, and a sub- 
sequent microwave pulse applied to the resonator is used to 
perform a joint dispersive read-out. 



used for joint three-qubit read-out by measuring its trans- 
mission [IS]. At the same time, it serves as a coupling 
bus for the qubits [H] . The qubits have a ladder-type en- 
ergy level structure with sufficient anharmonicity to allow 
individual microwave addressing of different transitions. 
We use the first two energy levels as the computational 
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qubit states, |0) and |1), and use the second excited state, 
1 2), to perform two-qubit and qubit-qutrit operations (a 
qutrit is a quantum ternary digit). From spectroscopy, 
we deduce a bare resonator frequency v T = 8.625 GHz 
with a quality factor of 3300; maximum qubit transi- 



tion frequencies v 
and ^c ax 



£ ax = 6.714 GHz, 



ix = 6.050 GHz 
4.999 GHz; and respective charging energies 
E c /h = 0.264, 0.296 and 0.307 GHz (h, Planck's con- 
stant) and qubit-resonator coupling strengths <?/27r = 
0.36, 0.30 and 0.34 GHz for qubits A, B and C. At the 
maximum transition frequencies, we find respective qubit 
energy relaxation times of T\ = 0.55, 0.70 and 1.10 lis 
and phase coherence times of T| = 0.45, 0.6 and 0.65 its 
for qubits A, B and C. 

In the conventional realization of the Toffoli gate, a 
not operation is applied to the target qubit (C) if the 
control qubits (A, B) are in the state |11). In our set-up 
it is more natural to construct a variation of the Toffoli 
gate shown in Fig. la in which the state of the target 
qubit is inverted if the control qubits are in |01). This 
gate can easily be transformed to the conventional Toffoli 
gate by a redefinition of the computational basis states 
of qubit A or by applying two 7r-pulses on qubit A. 

The Toffoli gate can be constructed from a 'controlled- 
controlled-phase' (CCPHASE) sandwiched between two 
Hadamard gates acting on the target qubit as shown in 
Fig. lb. A CCPHASE gate leads to a phase shift of 7r 
for state |1) of the target qubit if and only if the control 
qubits are in state 1 01) . In other words, this corresponds 
to a sign change of only one of the eight computational 
three-qubit basis states: 1 01 1) -H> — 1 01 1) . 

The basic idea of 'hiding' states by tranforming them 
into non-computational states to simplify the implemen- 
tation of a Toffoli gate was theoretically proposed in 
refs. l20l and l21l and has been experimentally implemented 
for linear optics and ion trap systems [5] [Uj. The im- 
plementation of the scheme of ref. [20 in our set-up 
would require three controlled-phase (CPHASE) gates, 
six single-qubit operations and two single-qutrit opera- 
tions. Instead, we construct the CCPHASE from a sin- 
gle two-qubit CPHASE gate and two qubit-qutrit gates. 
The latter gates are called tt-SWAP and 3tt-SWAP, re- 
spectively (Fig. lc, red frames). The application of a 



TABLE I. List of states after each step of the CCPHASE 
gate. The state 1 011) acquires a phase shift of 7r during the 
CPHASE pulse; the states 1 11a;) are transferred to i|20a;}, 
'hiding' them from the CPHASE gate; and the initial states 
\xQy) and |010) do not change during the sequence. 
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single CPHASE gate to qubits B and C (Fig. lc, blue 
frame) inverts the sign of both 1 1 1 1) and 1 1 1) . To cre- 
ate the CCPHASE operation, the computational basis 
state | 111) is transferred to the non-computational state 
i|201) by the 7T-SWAP gate, effectively hiding it from 
the CPHASE operation acting on qubits B and C. After 
the CPHASE operation, 1 1 1 1) is recovered from the non- 
computational level i|201) by the 37T-SWAP gate. Al- 
ternative approaches using optimal control of individual 
qubits for implementing a Toffoli gate in a single step 
have been proposed [55] and recently analyzed in the 
context of the circuit quantum electrodynamics architec- 
ture [23] . 

All three-qubit basis states show three distinct evolu- 
tion paths during our CCPHASE gate (Table 1). Only 
input state |011) is affected by the CPHASE gate acting 
on qubits B and C, which transfers |011) to the desired 
state, — 1 011) . The states |llx) with x € {0, 1} are trans- 
ferred by the tt-CPHASE gate to the states i\2Qx). The 
subsequent CPHASE gate therefore has no influence on 
the state. The last gate (3tt-CPHASE) transfers i\2Qx) 
back to | lis). Together the two SWAP gates realize a 
rotation by 4-7T, such that the state |llx) does not acquire 
any extra phase relative to the other states. The states 
of the last group (|010) and \x0y) with y € {0, 1}) do not 
change during the CPHASE gate sequence. 

The actual experimental implementation of the Tof- 
foli gate consists of a sequence of microwave and flux 
pulses applied to the qubit local control lines (Fig. Id). 
The arbitrary rotations about the x and y axis |24j are 
realized with resonant microwave pulses applied to the 
open transmission line at each qubit. We use 8-ns-long, 
Gaussian-shaped DRAG-pulses [24j [25] to prevent pop- 
ulation of the third level and phase errors during the 
single-qubit operations. Few-nanosecond-long current 
pulses passing through the transmission lines next to the 
superconducting loops of the respective qubits control 
the qubit transition frequency realizing z-axis rotations. 
All two-qubit or qubit-qutrit gates are implemented by 
tuning a qutrit non-adiabatically to the avoided cross- 
ing between the states 1 11a;) and |20x) or, respectively, 
|xll) and |a;20) (refs [PJJ [H and [37]). During this time, 
the system oscillates between these pairs of states with 
respective frequencies 2</-j^ and 2J^^ . With inter- 



action times 7r/(2J{^ ) 



7 ns, 3tt/(2J^ ) = 21 ns 
and 7r/(2J]^ ) = 23 ns, we realize a 7T-SWAP and a 
3tt-SWAP between qubits A and B and a CPHASE gate 
between qubits B and C, respectively. Our use of qubit- 
qutrit instead of single-qutrit operations allows for a more 
efficient construction of the Toffoli gate. Direct realiza- 
tion of the scheme proposed in ref. !20l in our system would 
require eight additional microwave pulses (used to imple- 
ment six single-qubit and two single-qutrit gates) with a 
twofold increase in overall duration of the pulse sequence 
with respect to our scheme. 
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FIG. 2. Truth table of the Toffoli gate. The state of qubit 
C is inverted if qubits A and B are in the state |01). The 
fidelity of the truth table is F = (l/8)Tr [U cxp U idc!i i} = 76.0% 



We have characterized the performance of this realiza- 
tion of a Toffoli gate by measuring the truth table, by 
full process tomography [28] and by Monte Carlo process 
certification [29l [30]. The truth table depicted (Fig. 2) 
shows the population of all computational basis states 
after applying the Toffoli gate to each of the computa- 
tional basis states. It reveals the characteristic properties 
of the Toffoli gate, namely that a NOT operation is ap- 
plied on the target qubit (C) if the control qubits (A 
and B) are in the state |0T). The fidelities of the out- 
put states show a significant dependence on qubit life- 
time. In particular, input states with qubit A (with the 
shortest lifetime) in the excited state generally have the 
worst fidelity, indicating that the protocol is mainly lim- 
ited by the qubit lifetime. The fidelity of the measured 
truth table, U cxp , with respect to the ideal one, t/jdeab 
namely F = (l/8)Tr [U c ^ p U idca i] = 76.0%, shows the av- 
erage performance of our gate when acting onto the eight 
basis states. 

As an essential addition to the classical characteriza- 
tion of the gate by the truth table, we have performed 
full, three-qubit process tomography and reconstructed 
the process matrix, Xoxp, to characterize the quantum 
features of the Toffoli gate completely, overcoming the 
limited characterization provided by measurements of the 
phase fidelity only [IB]. For this purpose, we prepared a 
complete set of 64 distinct input states by applying all 
combinations of single-qubit operations chosen from the 
set {id, 7r/2 x , n/2 y , ir x } for each qubit, and performed 
state tomography on the respective output states. The 
process matrix reconstructed directly from the data has 
a fidelity of F = Tr[xcx P Xidcai] = 70 ± 3% (the error 
represents a 90% confidence interval), where Xideai is the 



ideal process matrix. Using a maximum-likelihood pro- 
cedure [31] to correct for unphysical properties of Xcxp, 
we find that the obtained process matrix, x^xp: has a fi- 
delity of F = Tr[x^pXidcai] = 69 % with expected errors 
at the level of 3 %. In Fig. 3a, Xcxp shows the same key 
features as Xideai (Fig- 3b). 

To gain an accurate alternative estimate of the process 
fidelity without resorting to a maximum-likelihood proce- 
dure, we implemented Monte Carlo process certification 
following the steps described in ref. [29] First we define a 
Pauli observable as P n — Ylf=i gPnji a product of six 
single-qubit operators chosen from the set of the identity 
and the Pauli operators (p n ,j G {1, &x, &y, <?>})■ Then we 
determine the 232 observables with non-vanishing expec- 
tation values P n — TrlpxPn] 7^ 0, where pr is the Choi 
matrix of the Toffoli process. For each P n we prepare all 
(2 3 = 8) eigenstates of the product of the first three oper- 
ators comprising P n , apply the Toffoli operation to these 
states and measure the expectation value of the prod- 
uct of the last three operators in P n . Averaging over 
the results obtained with all eigenstates provides an esti- 
mate of P n . Extracting all 232 expectation values in this 
way allows us to estimate the fidelity of the Toffoli gate 
as 68.5 ± 0.5% using Monte Carlo process certification, 
which is in good agreement with the fidelity evaluated 
using tomography. 

The scheme that we use to implement the Toffoli gate 
is generic and can readily be applied to other systems be- 
cause the majority of the quantum systems used as qubits 
have additional energy levels at their disposal. Reduction 
of the total gate time by use of qubit-qutrit gates together 
with the recent advances in the extension of the coherence 
times of the superconducting circuits [32] [33J indicates a 
path towards the realization of practical quantum error 
correction. 

We thank Stefan Filipp, Alexandre Blais for useful 
discussions and Kiryl Pakrouski for his contributions 
in early stages of the experimental work. This work 
was supported by the Swiss National Science Foundation 
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